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ABSTRACT 

We attempt to explain the properties of the Becklin-Neugebauer (BN) object as a runaway B star, as 
originally proposed by Plambeck et al. (1995). This is one of the best-studied bright infrared sources, 
located in the Orion Nebula Cluster — an important testing ground for massive star formation theories. 
From radio observations of BN's proper motion, we trace its trajectory back to Trapezium star 9 1 C, the 
most massive (45M ) in the cluster and a relatively tight (17 AU) visual binary with a B star secondary. 
This origin would be the most recent known runaway B star ejection event, occurring only ~ 4000 yr 
ago and providing a unique test of models of ejection from multiple systems of massive stars. Although 
highly obscured, we can constrain BN's mass (~ 7Mq) from both its bolometric luminosity and the 
recoil of 6 C. Interaction of a runaway B star with dense ambient gas should produce a compact wind 
bow shock. We suggest that X-ray emission from this shocked gas may have been seen by Chandra: 
the offset from the radio position is ~ 300 AU in the direction of BN's motion. Given this model, we 
constrain the ambient density, wind mass-loss rate and wind velocity. BN made closest approach to the 
massive protostar, source "I" , 500 yr ago. This may have triggered enhanced accretion and thus outflow, 
consistent with previous interpretations of the outflow being a recent (~ 10 3 yr) "explosive" event. 

Subject headings: stars: formation — stars: kinematics — stars: winds, outflows 



1. INTRODUCTION 

The Becklin-Neugebauer (BN) object was discovered as 
a bright 2 fj,m source (Becklin & Neugebauer 1967), about 
0.1 pc (45") in projection from the Trapezium stars of the 
Orion Nebula Cluster (ONC) (e.g. O'Dell 2001), about 
450 pc distant (Genzel & Stutzki 1989). It was initially 
thought to be a protostar (i.e. embedded and in main 
accretion phase), since it is infrared bright but optically 
faint. The extinction-corrected spectrum peaks at ~ 5/im 
(T ~ 550 K), with Ay ~ 17 mag (from depth of sili- 
cate 9.8 /mi absorption), yielding a bolometric luminosity 
~ 25OOZ/0 from a "core" region (~ 1") and ~ 1O 4 L from 
an extended region (~ 10") (Gezari, Backman, & Werner 
1998). These estimates are uncertain as they are derived 
from single temperature blackbody fits to limited portions 
of the spectrum: as Gezari et al. describe, the BN re- 
gion probably exhibits a range of temperatures. The ex- 
tinction correction is also uncertain: the above estimate 
agrees with Scoville et al.'s (1983) determination from the 
Bra/Br7 emission line ratio assuming case B recombina- 
tion theory. However Bunn, Hoare, & Drew (1995) argued 
from the Pf7 line strength that case B theory is not appli- 
cable and Ay > 18 mag by a significant amount. 

If the 2500 — 10 4 L Q corresponds to the luminosity of a 
solar metallicity zero age main sequence star, then this has 
mass ~ 8.4 — 12M , radius r* ~ 3.4 — 4.3 Rq and sur- 
face temperature ~ 22, 000 - 28, 000 K (Schaller et al. 
1992). BN has 8.3 day periodic variability with ~ 0.2 mag 
peak to trough change at H and K (Hillenbrand, Carpen- 
ter & Skrutskie 2001). This behavior is similar to that 
of JW 660, a mid-B spectral type (~ 6M ) star also in 
the ONC with a 6.15 day period (Mandel & Herbst 1991). 
These periods are somewhat longer than those typical of 
slowly pulsating B stars (Mathias et al. 2001). 



BN is detected in the radio (Moran et al. 1983; Garay, 
Moran & Reid 1987; Churchwell et al. 1987; Menten & 
Reid 1995; Plambeck et al. 1995) and just resolved at 2 cm 
with diameter 31 x 18 AU (0.07" x 0.04") (Churchwell et al. 
1987). The spectrum rises as i/ 125 to ~ 100 GHz and then 
may flatten. Plambeck et al. (1995) present a non-unique 
model of this emission: a 10 4 K H II region with uniform 
density core (n e = 5 x 10 7 cm~ 3 ; radius r = 7 AU) and 
steeply declining envelope with n e cx (r/ro 
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2. THE MOTION OF BN 

From radio data with a time baseline of 8.7 years (1986.3 
- 1995.0), Plambeck et al. (1995) reported a proper mo- 
tion of BN relative to source "I" (a massive protostar 
forming from the Orion hot core, Greenhill et al. 1998; 
Tan 2003) of - 0.02" yr" 1 (~ 50 km s" 1 in the sky plane 
for BN at 450 pc, adopted throughout), towards the NW. 
From the four data points of Plambeck et. al., we find 
0.021 ± 0.004" y^ 1 (45 ±8 km s _1 ) towards position angle 
(P.A.) —33 ± 11°, where errors are la assuming the indi- 
vidual coordinates are known to ~ 0.02". The relative 
positions have continued to be determined with BIMA 
(Berkeley-Illinois-Maryland-Array) (Plambeck 2004, pri- 
vate comm.). From a joint analysis of the full data 
set (extending to 2004) we estimate a proper motion of 
0.0181 ± 0.0022" yr" 1 (38.7 ± 4.7km s" 1 ) towards P.A. 
—37.7 ± 5°. Following the velocity vector backwards as- 
suming no acceleration, we see BN was near the Trapez- 
ium stars ~ 4000 yr ago (Fig. 1) and made a close (< 2") 
projected passage of source "I" ~ 500 yr ago. 

Scoville et al. (1983) observed several spectroscopic 
tracers of BN and its circumstellar nebula, reporting a ve- 
locity relative to the local standard of rest of +21 km s _ . 
This is significantly larger than the ONC mean (+8 ± 
3.5 km , Walker 1983) and the gas of the Orion hot 
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core and larger scale molecular cloud (+9 km s~ , Genzel 
& Stutzki 1989). Thus relative to the Orion cloud and star 
cluster, BN is moving at «bn — 41 ± 6 km s _1 . 

If BN originated near the Trapezium center then it is 
now ~ 0.05 pc further away from us. BN's extinction 
(A v > 17, Gezari et al. 1998, Bunn ct al. 1995), im- 
plies a large intervening column of gas and dust, Ah > 
6 x 10 22 cm~ 2 , 90% of which is behind the Trapezium stars. 
Wen & O'Dell (1995) considered the creation of an H II 
region by 9 1 C and its interaction with the nearby neutral 
gas. They argue this gas is ~ 0.15 pc behind the Trapez- 
ium stars for a uniform density model and somewhat closer 
for exponential models. The possible discrepancy of these 
estimates with the ^0.05 pc inferred from BN's kinemat- 
ics warrants further study. For a neutral layer thickness 
of < 0.05 pc, the mean density on the BN sight line is 
> 4 x 10 5 cm -3 . This material is likely associated with 
the Orion hot core that harbors source "I" . This core may 
have had an initial mass ~ 6OM (Tan 2003) and density 
at its surface n H , s ^ 10 6 cm~ 3 (McKee & Tan 2003), in 
reasonable agreement with the above estimate. 

3. CREATION OF A RUNAWAY: EJECTION FROM 6 1 C 

Two models are thought to explain runaway OB stars. 
Zwicky (1957) and Blaauw (1961) proposed that they orig- 
inate when a supernova occurs in a close binary. Alterna- 
tively, Poveda, Ruiz, & Allen (1967) hypothesized origin 
by dynamical ejection of a star from a multiple system. 
Hoogerwerf, de Bruijne & de Zeeuw (2001) considered past 
trajectories of a large sample (~ 50) of nearby (< 1 kpc) 
runaways and several pulsars. They found one case of a 
pulsar and B star that supports the supernova scenario. 
They also identified two B stars and a massive eccentric 
binary with an ejection event from the ONC ~ 2.5 Myr 
ago, which supports dynamical ejection (see also Gualan- 
dris, Portegies-Zwart & Eggleton 2004). 1 

BN, if created in the ONC, must represent an example 
of dynamical ejection, since the cluster is young (< 3 Myr; 
Palla & Stahler 1999) and there is no evidence for a recent 
supernova. Dynamical modeling (e.g. Leonard & Duncan 
1990) suggests the runaway should often leave behind a 
hard, eccentric binary of much greater total mass and an 
escape velocity from a location typical of the secondary's 
orbit that is comparable to the ejection velocity. 

All the Trapezium stars lie close to the possible past 
trajectories of BN (Fig. 1). Their binary properties are 
(Preibisch et al. 1999; Schertl et al. 2003): 9 1 A (B0, 
16M ) has a visual companion at 100 AU (4M Q ) and a 
spectroscopic companion at ~ 1 AU (~ 3M Q ); 9 1 B (7M Q ) 
has a spectroscopic companion at < 0.1 AU (<~ 2M Q ) and 
3 other more distant companions of at most a few solar 
masses; 9 1 C (06, 45M Q ) has a visual secondary at 17 AU 
(> 6M Q ), perhaps on an eccentric orbit (Schertl ct al. 
2003) and able to explain evidence for a spectroscopic com- 
panion (Donati et al. 2002; Vitrichenko 2002); 9 1 D (B0. 5) 
has no detected companions. We note 9 2 A (09.5, 25M ), 
the second most massive ONC star and located at (-130", 
-150") relative to source "I" , also lies along an extrapola- 
tion of BN's past trajectory. It has a spectroscopic com- 

1 This event implies that massive star formation has been occurring in the ONC for ~ 2.5 Myr, and that a large fraction (almost 50%) of 
massive (m* > IOMq) stars have been ejected. As Hoogerwerf ct al. (2001) point out, this affects determinations of the high-mass initial mass 
function from studies of star clusters. It also has implications for models of star cluster formation (e.g. Tan & McKee 2001), which assess the 
importance of feedback from massive stars: this is reduced if such stars are efficiently ejected from clusters. 



panion at 0.47 AU (~ 9M Q , e = 0.33) and a visual com- 
panion at 173 AU (7M Q ) (Preibisch et al 1999; Abt, Wang, 
& Cardona 1991). The escape velocity from the 9 1 C sys- 
tem from a location characteristic of the secondary's or- 
bit is 67(m i c /5OM ) 1 / 2 (r/2OAU)- 1 / 2 km s _1 , while it is 
~ 350 km s _1 from the inner 9 2 A system. These appear 
to be the only ONC stars capable of ejecting BN. 

If no other star was ejected with BN, then 
the parent system should be recoiling with a 
proper motion, evaluated for the case of 9 1 C, of 
3.6mas yr~ 1 (mBN/lOM0)(mgi C 7/5OM0)~ 1 , where lmas = 
0.001". This is several times larger than the dispersion of 
bright ONC stars (±0.70 ± 0.06 mas/yr, van Altena et al. 
1988), and much larger than the velocity expected if 9 1 C 
were in equipartition with the other cluster stars. In fact 
van Altena ct al. (1988) report a proper motion of 9 1 C 
of 2.3 ± 0.2 mas/yr towards PA. +142.4 ± 4° (corrected 
to J2000), with the quoted errors resulting from assum- 
ing ±0.15 mas/yr in the x and y velocity vectors (van 
Altena et al. 1988). The uncertainties in 9 1 C's motion 
are ~ 20% larger due to its brightness (van Altena 2004, 
private comm.), but this does not significantly change our 
results. 9 2 A is moving towards —30° at ~ 2.7 mas/yr, i.e. 
roughly towards the Trapezium and BN, so it is unlikely to 
have ejected BN. Note that these motions are relative to 
a frame defined by the mean proper motion of van Altena 
et al.'s sample of ONC stars, while BN's motion has been 
measured relative to source "I" . We are assuming that "I" 
and the ONC have negligible relative proper motion. The 
direction of 9 1 C , s motion is consistent with being exactly 
opposite to BN's (Fig. 1). We searched van Altena et al.'s 
sample for other relatively high proper motion stars that 
may have been involved in the ejection event: none were 
found. Assuming BN is the only ejected star and the initial 
system had negligible velocity with respect to the ONC, 
then BN's mass is tobn = 6.35 ± l.Q(mgic /§QMq)Mq. 
A ±0.7 mas/yr uncertainty in the initial motion of the 
system translates into an additional ±2M Q uncertainty in 
?tibn- BN is now ~ 10" away from source "I", causing a 
velocity gain of ~ 4.0(m//40M Q ) 1 / 2 km s _1 . Accounting 
for this boosts BN's mass estimate by ~ 10%. Treating BN 
as a massless test particle, its deflection angle due to "I" 
is 1.56°(TO/,»/20M Q )(6/1000AU)- 1 (w B N/36km s" 1 )- 2 , 
where b is the initial impact parameter and we now nor- 
malize to the expected protostellar mass of "I" (Tan 2003). 

BN passed close to source "I" about 500 yr ago. Our fa- 
vored interpretation is that BN did not form from the gas 
associated with "I" and the hot core (or Kleinmann-Low 
nebula) . This possibility was first suggested by Zuckerman 
(1973). However, we cannot rule out ejection from "I". A 
close passage of BN with source "I" and its protostellar 
disk may have led to enhanced angular momentum trans- 
port and accretion via tidal torques (e.g. Ostriker 1994). 
The protostellar outflow rate is predicted to be propor- 
tional to the accretion rate in a range of theoretical mod- 
els (e.g. Shu et al. 2000; Konigl & Pudritz 2000). Allen 
& Burton (1993) have argued that the outflow from "I" 
(i.e. from the Kleinmann-Low nebula) shows characteris- 
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tics consistent with an "explosive" , i.e. impulsive, event 
occurring about one flow crossing time ago (~ 10 3 yr). 

For a 3-body interaction, the kinetic energy of BN and 
6 1 C (9.1xl0 4 ^TOBN/7M Q )(w B N/36km s" 1 ) 2 crgs + 1.24 x 
10 46 (m e i c /50M Q )(u e i c /5km s" 1 ) 2 ergs) should not be 
significantly greater than the binding energy of the binary 
(1.16 x lO 47 (mi/45M Q )(m 2 /5M ) ergs, assuming e = 0). 
This comparison suggests that either the binary is quite 
eccentric (e.g. Schertl et al. 2003) or 9 1 C's component 
masses are larger than the adopted values. The period of 
the 9 1 C binary is about 10 years, so about 400 orbits have 
occurred since the ejection. Given a sufficiently accurate 
determination of BN and 6 1 C's motion and the 1 C binary 
orbit, details of the ejection could be constrained. 

4. THE STELLAR WIND BOW SHOCK 

A star emitting a wind with mass-loss rate m w and ter- 
minal speed v w and moving through a medium of density 
p a with velocity v* produces a bow shock at distance r^ s 
in the direction of motion. Assuming no penetration of 
ambient material into the wind bubble and balancing ram 
pressures, PwV^ = m w v w / (Airr^) = p a v^, we have 




assuming a hydrogen mass fraction of 0.7. 

X-ray emission from BN is offset by < 1" (Garmire 
et al. 2000; Feigelson et al. 2002; Feigelson and the 
Chandra Orion Ultradeep Project [COUP] 2004, private 
comm.). The most recent determination from COUP is 
offset from the 2003 radio position by 0.60 ± 0.1" towards 
P.A. -32 ± 14° (Fig. 2). The density of X-ray sources in 
this region is <~ 0.02 arcsec~2, so the probability of finding 
an unrelated source within l"of BN is ~ 0.06. The align- 
ment of the observed source with BN's direction of motion 
suggests that it may be produced at the wind bow shock. 
If the X-rays come from the inner edge of the thin shell 
of shocked wind and ambient gas, with each solid angle of 
wind material from a region about l"across contributing 
equally, then we can use the solution of Wilkin (1996) for 
the shell geometry to estimate rb s - The offset of the center 
of the X-ray emission is 90% of rb s ■ For the fiducial geom- 
etry, projection effects increase it by ~ 5%, so we estimate 
rbs = 300 ±50 AU. We note that in this interpretation the 
bow shock is much further from the star than the extent of 
the ionized region, which is tens of AU. The stellar ioniz- 
ing flux is confined within the wind even before it reaches 
the bow shock (c.f. Van Buren et al. 1990). 

The X-ray spectrum constrains the initial wind veloc- 
ity, since the temperature of the postshock (7 = 5/3) gas 
is kT = (3/16)^ = 1.96(^/m H )(Vl0 3 km s" 1 ) 2 keV, 
and the shock speed v s ~ v w . Based on 45 counts, 
Feigelson et al. (2002) reported that the emission is hard 
(kT > 10 keV) and heavily absorbed (7V H ~ 4x 10 22 c m - 2 ). 
The observed (0.5-8 keV) luminosity is 2.5 x 10 29 crgs s _1 , 
becoming ~ 4 x 10 29 ergss^ 1 once corrected for absorp- 
tion. With so few counts the uncertainties in these pa- 
rameters are large: ~ 50% in kT and a factor of ~3 in 
iVn- Nevertheless these results suggest the presence of a 



fast wind. The escape speed from main sequence stars 
with to* = 8.4 — 12M (described in §1) are v csc ^ = 
970 — 1030 km s _1 . Consistency with the observed X-ray 
spectrum may be possible for a model with a somewhat 
lower temperature (i.e. ^few keV), an increased amount 
of absorption, and a much higher luminosity (see below). 

The wind bow shock model predicts a relatively con- 
stant X-ray source, unless there are large variations in the 
power of the stellar wind or the density of the ambient 
medium. The limited variability reported by Feigelson et 
al. (2002) is consistent with a constant source. The model 
also predicts that the source may show extended emission 
on scales ~ 1" but perhaps of low surface brightness. 

The mass-loss rates from young B stars are uncer- 
tain. Even for stars on the main sequence the theoretical 
mass-loss rate due to a line-driven wind is uncertain be- 
cause the star may be close to the "bi-stability jump" at 
T* ~ 21000 K, where there is a change in the ionization 
state of the lower wind layers near the sonic point (Vink 
et al. 1999). Using the results of Vink et al. (2001), the 
m* = 8.4M case can have m w ~ 2 — 70 x 10 _11 M Q yr -1 , 
with v w ~ 2.7 — 1.3f C sc,*- F° r t ne m * = 12M Q case 
m w = 1.15 x 10~ 9 M Q yr -1 and v w ~ 1.3u esc ,*- 

It takes a time £kh = a ff /3GTO 2 /(2r*L*) = 1.57 x 
10\ P(m.J IOMq^^/IORq)- 1 ^*/ 10* Lq)- 1 yr to set- 
tle to the main sequence, where (3 is the mean ratio of gas 
to total pressure (~ 1 for to* ~ 10M Q ) and a g = 3/(5 — n) 
for polytropcs with n < 5. Approximating n = 3 then 
£kh — 20— 8x 10 4 yr for the to* = 8.4— 12M Q cases, respec- 
tively. Thus it is possible that the star is still in a pre-main 
sequence phase. If so, then empirically we expect higher 
mass-loss rates, ~ 1O _8 -1O _6 M yr^ 1 (Nisini et al. 1995) 
and somewhat lower wind velocities (by factors of a few), 
compared to main sequence wind models. Indeed mod- 
els of slower, denser winds have been proposed to explain 
radio continuum emission and near-infrared hydrogen re- 
combination lines observed from BN. Scoville et al. (1983) 
estimated m w ~4x 10~ 7 M Q yr -1 and v w > 100 km s _1 . 
Hoflich & Wehrse (1987) presented a non-LTE model with 
m w = 3 x 10 _7 M Q yr" 1 and v — 20km s~ . Bunn 
et al. (1995) observed Br 7 line wings extending to 
~ 220 km s -1 . Note that the wind luminosity is L w = 
3.15x1O 34 (to„,/1O- 7 M yr- 1 )(^/1000km s~ 4 ) 2 ergs s" 1 
(= 8.2L Q ), which is small compared to the bolometric lu- 
minosity, but large compared to the X-ray luminosity. 

Given the uncertainties in fn w and v w , we cannot use the 
bow shock model to estimate unambiguously the ambient 
density. However, the product of these three quantities is 
constrained (eq. 1). Estimates of the mean density (e.g. 
§2) are about an order of magnitude higher that the fidu- 
cial value shown in eq. 1 , suggesting (in the context of our 
bow shock model) that rh w and/or v w are relatively large 
compared to the adopted fiducial values. 

We note that some of the H recombination line emission 
that motivated "slow" wind models may be generated by 
shocked ambient interstellar material. From mass continu- 
ity of gas entering the shocked shell and flowing around the 
wind bubble, the shell's column density is approximately 
n H , a r bs = 4.5 x 10 20 (n H , a /10 5 cm~ 3 )(r bs /300 AU) cm" 2 . 
As an aside, this is small compared to the total inferred 
absorbing column, which must be dominated by inter- 
vening ambient material. For an isothermal shock with 
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T a = 100 K and fj, a = 2.35mn so that c a = 0.59 km s~ , 
then the Mach number is ~ 70 and the compression ratio 
~ 5000 for v* = 40 km s _1 . Thus an upper limit to the 
shocked density is ~ 10 10 cm -3 , assuming nu. 

a ~ 10 cm 

The actual densities may be much smaller if magnetic 
fields contribute significantly to the total pressure. 

From CO emission features, Scovillc et al. (1983) 
reported the presence of molecular gas with jih 2 ~ 
10 7 - 10 12 cm" 3 and T ~ 600 - 3500 K. Some of this 
emission may be due to the swept-up shell, but the highest 
densities are difficult to explain with the bow shock model. 
A more likely alternative is that the densest gas, < 1 AU 
from the star, is the remnant of an accretion disk. This is 
the favored interpretation to explain CO emission from a 
number of luminous young stellar objects (e.g. Ishii et al. 



2001). The disk is probably not being fed significantly by 
swept-up material: the Bondi-Hoyle accretion rate, ignor- 
ing stellar feedback, is still relatively small, tobh = 1-3 x 
10- 9 (n H: a/10 5 cm- 3 )(m„/10M Q ) 2 (w st /40km s _1 )- 3 M yr" 1 
(e.g. compared to the expected mass loss rate). If, how- 
ever, the ambient medium is clumped on very small scales 
then it would be easier to occasionally bring dense inter- 
stellar molecular gas close to BN at higher rates. 
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Fig. 1. — The velocity vector and inferred past trajectory of BN (central dotted line), based on radio observations (Churchwell et al. 1987; 
Menten & Reid 1995; Plambeck et al. 1995; Plambcck 2004, private comm.). Coordinates are relative to source "I". An approximate estimate 
of the la uncertainty in the position angle of the trajectory is shown by the outer dotted lines. BN was near the Trapezium stars, about 
4000 yr ago. In particular, wc argue (see text) that it was ejected from the 8 1 C binary system, the proper motion of which is also shown. 
Note that BN made closest passage 500 yr ago to source "I" , a massive protostar. The dashed line shows the approximate size of an initial 
equilibrium 60Mq core (McKee & Tan 2003; Tan 2003), bounded by a (nonthermal) pressure ~ GS 2 with E = 1 g cm -2 . 
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Fig. 2. — Offset of X-ray emission (Feigelson et al. 2004, private comm.) from BN (Jan 2003, Plambeck, private comm.). The coordinates 
are relative to source "I". The size of the crosses indicate positional uncertainties. Dotted lines show past and future trajectories of BN, as in 
Fig. 1 (normalized to pass through the mean position of the data from 1986.3 to 2004: hence slight offset from 2003 position). The dashed 
line shows the projected extent of the bow shock solution of Wilkin (1996) for a minimum offset of r bs = 300 AU (0.67"). 



